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Abstract 



Let B be a point robot moving in the plane, whose path is constrained to forward motions with 
curvature at most one, and let P be a convex polygon with n vertices. Given a starting configuration 
(a location and a direction of travel) for B inside P, we characterize the region of all points of P 
that can be reached by B, and show that it has complexity 0(n). We give an 0(n 2 ) time algorithm 
to compute this region. We show that a point is reachable only if it can be reached by a path of type 
CCSCS, where C denotes a unit circle arc and S denotes a line segment. 

1 Introduction 

The problem of planning the motion of a robot subject to non-holonomic constraints [16, 25 (for instance, 
bounds on velocity or acceleration [5] HU [22] > bounds on the turning angle) has received considerable 
attention in the robotics literature. Theoretical studies of non-holonomic motion planning are far sparser. 

In this paper we consider a point robot in the plane whose turning radius is constrained to be at 
least one and that is not allowed to make reversals. This restriction corresponds naturally to constraints 
imposed by the steering mechanism found in car-like robots. We assume that the robot is located at a 
given position (and orientation) inside a convex polygon, and we are interested in the set of points in 
the polygon that can be reached by the robot. We put no restriction on the orientation with which the 
robot can reach a point. 

The lack of such a restriction distinguishes our work from most of the previous theoretical work on 
curvature-constrained paths, which usually assumes that not a point, but a configuration (a location with 
orientation) is given. Dubins [12] was perhaps the first to study curvature-constrained shortest paths. 
He proved that a curvature-constrained shortest path from a given starting configuration to a given final 
configuration consists of at most three segments, each of which is either a straight line or an arc of a 
unit-radius circle. Reeds and Shepp |21j extended this characterization to robots that are allowed to 
make reversals. Using ideas from control theory, Boissonnat et al. [B] gave an alternative proof for both 
cases, and Sussmann [55] extended the characterization to the 3-dimensional case. 

In the presence of obstacles, Fortune and Wilfong [T3] gave a single-exponential decision procedure to 
verify if two given configurations can be joined by a curvature-constrained path avoiding the polygonal 
obstacles. On the other hand, computing a shortest bounded-curvature path among polygonal obstacles 
is NP-hard, as shown by Reif and Wang [53]. Wilfong [37] designed an exact algorithm for the case 
where the curvature-constrained path is limited to some fixed straight "lanes" and circular arc turns 
between the lanes. Agarwal et al. [2] considered the case of disjoint convex obstacles whose curvature 
is also bounded by one, and gave efficient approximation algorithms. Boissonnat and Lazard [5] gave 
a polynomial-time algorithm for computing the exact shortest paths for the case when the edges of the 
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obstacles are circular arcs of unit radius and straight line segments. Boissonnat et al. [7] gave an 0(n 4 ) 
algorithm for finding a convex and simple path of bounded curvature within a simple polygon. Agarwal et 
al. [1] presented an 0(n 2 logn)-time algorithm to compute a curvature-constrained shortest path between 
two given configurations inside a convex polygon. They also showed that there exists an optimal path 
that consists of at most eight line segments or circular arcs. For general polygonal obstacles, Backer and 
Kirkpatrick [5] recently gave the hrst complete approximation algorithm, improving on earlier work that 
approximated the shortest "robust" path [T3 [3] . 

At least two interesting problems have been studied where not configurations but only locations 
for the robot are given. The first problem considers a sequence of points in the plane, and asks for 
the shortest curvature-constrained path that visits the points in this sequence. In the second problem, 
the Dubins traveling salesman problem, the input is a set of points in the plane, and asks to find a 
shortest curvature-constrained path visiting all points. Both problems have been studied by researchers 
in the robotics community, giving heuristics and experimental results [2H HH1 [H] . From a theoretical 
perspective, Lee et al. |17j gave a linear-time, constant-factor approximation algorithm for the first 
problem. No approximation algorithms are known for the Dubins traveling salesman problem. 

Our result is a characterization of the region of points reachable by paths under curvature constraints 
from a given starting configuration inside a convex polygon P. We show that all points reachable from 
the starting configuration are also reachable by paths of type CCSCS, where C denotes an arc of a 
unit-radius circle and S denotes a line segment. When P has n vertices, we show that the reachable 
region has complexity 0(n), and we give an 0(n 2 ) time algorithm to compute this region. 

2 Terminology and some lemmas 

Let P be a convex polygon in the plane. A configuration s = (s, d) is a point s together with a direction 
of travel d (a unit vector). By a path, we mean a continuously diffcrcntiable curve (the image of a 
C^-mapping of [0, f] to K 2 ) with average curvature bounded by one in every positive-length interval. 
Unless stated otherwise, we assume that a path is completely contained in P. A configuration on the 
path 7r is a configuration s = (s, d) with s on 7r such that d is the forward tangent to 7r in s. The starting 
configuration of 7r is the starting point of 7r with its forward tangent. A simple path is a path with no 
self-intersection; we allow the endpoints of a simple path to coincide, in which case we call it a simple 
closed path. (Hence, a simple closed path is smooth except possibly at one point.) 

Given a configuration s = (s,d), the left disk -Dl(s) (right disk Dr(s)) is the unit disk touching s 
and completely contained in the left (right) halfplane defined by the directed line through d. (All unit 
disks in this paper are unit-radius disks.) 

The left directly accessible region lda(s) is the set of all points in P that can be reached by a path 
with starting configuration s consisting of a single (possible zero-length) circular arc on the boundary 
of Dl(s) followed by a single (possible zero-length) line segment. The right directly accessible region 
is defined analogously. The directly accessible region da(s) is the union of the left and right directly 
accessible region. 

Pestov-Ionin lemma. The following lemma is perhaps the foundation for all our results. In a slightly 
less general form, it was proven by Pestov and Ionin [20] . Recently, it has been used for a curve 
reconstruction problem [14]. 

Lemma 1 (Pestov-Ionin). Any simple closed path contains a unit disk in its interior. 

For sake of completeness, we sketch a proof analogous to Pestov and Ionin's. 

Lemma 2. Let D be a closed disk, and T be a simple path with endpoints (a, b) such that TD D = {a, b}. 
Then there is a unit disk touching T \ {a, b} that lies within the region 1Z bounded by T and the exterior 
arc of 3D (See Figure [7]). 

Proof. We proceed by induction on the length I of T (or, more precisely, by induction on L^/ttJ). 

When £ < n, we can prove by integration that a unit disk tangent to T does not cross T. We consider 
a unit disk Do tangent to V at m ^ {a, b} on the interior side. Since Do \ D has only one connected 
component and m G Do \ D, clearly Do is contained in 1Z. 
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Figure 1: Illustration of Lemma [2] The region 1Z is shaded. 

Otherwise, let m be the point halfway between a and b on T. Let D\ be the largest disk contained 
in 1Z that is tangent to T in m. If the radius of Di is larger or equal to 1, then we are done. If not, Di 
must touch dlZ in another point m' . Clearly m! lies on F, and the length of the arc F' of T between m 
and m! is at most i/2. By the induction hypothesis, a unit disk D2 lies inside the region IZi bounded 
by r' and D\. Since 1Z\ C 1Z, the lemma follows. □ 

Lemma [l] follows from Lemma [2] by observing that it still holds when D degenerates to a point. 

Filling. By fil(P) we denote the set of all unit-radius disks that are completely contained in P, and 
we let Fil(P) be the union of all the disks in fil(P). (See Figure [2]) Both fil(P) and Fil(P) will be 
called the filling of P. 




Figure 2: The filling Fil(P) is shaded in light grey, and the three pockets are shaded in dark grey. 



Pockets. The connected components of P\Fil(P) are called the pockets of P. A pocket of P is bounded 
by a single circular arc (lying on one disk of fil(P)) and a connected part of the boundary of P. The 
first and last edge on this connected chain are called the mouth edges of the pocket. Its extremities are 
called the mouth points. The mouth edges form an angle smaller than tt (this is equivalent to observing 
that the mouth points lie on the same disk of fil(P) and form an angle smaller than ir) [T]. Agarwal et 
al. [1] proved the following lemma. 

Lemma 3 (Pocket lemma). A path entering a pocket from Fil(P) cannot leave the pocket anymore. 

Proof. We consider a pocket K bounded by a disk D. Suppose that there is a path T that enters and 
leaves K. There is a subpath V of V whose endpoints lie on D and whose other points are in the interior 
of K. By Lemma [2j there is a unit disk D' that touches T' at a point in the interior of K, and is 
contained in K U D. Hence, D' £ fil(P), and D' intersects the interior of K, a contradiction. □ 

Reachability for a union of disks. For a set T> of unit disks, we let CONv(2?) denote the set of all 
unit disks contained in the convex hull of 1J T>. Equivalently, CONv(2?) consists of the unit disks centered 
at points of the convex hull of the set of all centers of the disks in T>. 

Observation 4. Let V be a set of unit disks in the plane. Then p| T> = p| CONv(2?). 
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Given a convex set Q, a configuration on the boundary of Q is a configuration s = (s, d) with s on 
the boundary of Q and such that d is tangent to the boundary of Q in s. 

Lemma 5. Let T> be a set of unit disks, and let s be a starting configuration on the boundary of 
IJ CONv(2?) . Then no point in the interior of |~) T> can be reached by a path starting at s and contained 
in (J CONv(2?) ; or even in any convex polygon P such that fil(P) = CONv(2?). 

Proof. Assume to the contrary that there is a path 7 with starting configuration s on the boundary of 
(J CONv(X>) and ending point t in the interior of f] D. Assume for the moment that 7 lies completely in 

lJCONV(X>). 

We extend 7 to infinity using a straight ray, such that the extended path is still C 1 . We then extend 
7 backwards, by attaching a single loop around the boundary of (J CONV(P) at s. To summarize, the 
extended path 7' starts at s, makes a single loop around the boundary of (J CONv(D), then follows the 
original path 7, and finally escapes to infinity along a straight line. We can now construct a simple, closed 
path 7" as follows: starting at infinity, we follow 7' backwards, until we encounter the first intersection 
of 7' with the part that we have already seen. Such an intersection must exist since the two extensions 

^ \ S 
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Figure 3: Proof of Lemma [5j The set T> consists of three disks. The convex hull CONV(2?) is shaded in 
light grey, and f] T> is in dark grey. The path 7 from s to t is enlarged into the path 7" . 

intersect. We define 7" to be the part of 7' between these two self-intersection points. Observe that t 
lies either on or outside the closed loop 7". 

By the Pestov-Ionin Lemma, 7" contains a unit-disk D. Since 7" is contained in (Jconv(2?), we 
have D £ CONV(Z>). Consequently, the interior of f]T> lies in the interior of D, a contradiction with the 
fact that t must lie on or outside 7". 

The lemma still holds when we allow the path to lie inside a convex polygon P with fil(P) = 
CONV(P). After all, by Lemma |3j the path cannot return to Fil(P) after it has left it. □ 

The characterization. Consider a starting configuration s on the boundary of the filling Fil(P). It 
is easy to see that any point in P not in p| fil(P) can be reached by a path from s. On the other hand, 
by Lemma [5] no point in p|FlL(P) can be reached, and so we have a complete characterization of the 
region reachable from s as the complement of P|fil(P). 

If the starting configuration lies on the boundary of an arbitrary unit disk contained in P, the same 
characterization holds. For arbitrary starting configurations, however, the situation becomes far more 
complicated. There is the possibility that no path starting at s is tangent to the boundary of Fil(P), 
and the filling has no relation to the reachable region. 

If there exists a path starting at s that is tangent to the boundary, all points outside |~| fil(P) are 
reachable, but it is still possible that some points inside f] fil(P) are reachable, for instance because they 
lie in the directly reachable area da(s) = lda(s) U RDA(s), or because s lies in a pocket with additional 
maneuvering space (that we would not have been able to exploit if starting inside Fil(P) by Lemma [3| . 
(See Figure [I]). 

In the rest of this paper, we give a complete characterization of the reachable region, for any starting 
configuration in P. Let us denote the set of points t g P such that t is reachable by a path starting from 
a configuration s by reach(s). 
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Figure 4: Points in the grey region are reachable from s, but are neither in da(s) nor in the complement 

of D = f|FIL(P). 

3 Paths starting along the boundary 

In this section, we assume that the starting configuration s is on the boundary of P (recall that this 
means also that the direction is tangent to the boundary). Without loss of generality, we also assume 
that the direction of s is counterclockwise along the boundary, so that points of P are reached locally 
by a left turn from s. It turns out that in this situation we can restrict ourselves to paths containing no 
right-turning arcs. 

The forward chain fc(s) is the longest subchain of the boundary of P, that starts counterclockwise 
from s, and that turns by an angle at most tt. (See Figure pi) In other words, when s is directed 




Figure 5: The forward chain Fc(s). 

vertically upward, this chain contains all the edges of P that are above its interior, as well as the part of 
the edge that contains s and is above s, and, if there is one, the other vertical edge of P. 

If the forward chain intersects the interior of Z?l(s), then we have the following simple description of 
the reachable region. 

Lemma 6. If the forward chain Fc(s) intersects the interior of Dl(s), then REACH(s) = lda(s). 

Proof. The left directly accessible region lda(s) can be enlarged to a pocket of the left disk Dl(s). Thus, 
no point outside lda(s) is reachable by the Pocket Lemma (Lemma [3]). □ 

We denote by LFIl(s) the set of disks contained in P U Dl(s) that touch the forward chain. Note 
that lfil(s) always contains the left disk Dl(s). Let us remark that the set of centers of the disks in 
lfil(s) need not be connected. For example, Figure [6] shows a situation where lfil(s) consists of just 
three disks; their centers are marked by black dots (in general, any number of connected components is 
possible). For a disk F £ lfil(s), we define lda(P) as LDA(r), where r is the first configuration on 
fc(s) touching F (by the above, this is well defined). Note that if F £ fil(P), then lda(P) is simply 
the complement of the interior of F. We continue with a lemma on the reachable points outside lda(s). 

Lemma 7. Let P and s be as above, with D^(s) ^ fil(P), and the forward chain fc(s) not intersecting 
the interior of Dl(s). Suppose that there is a path 7 from s to t, where t (£ lda(s). Then there exists a 
disk F such that F £ fil(P) and t is not in the interior of F , or there exists a disk F £ lfil(s) such 
that t £ lda(P). 
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Figure 6: Example where the left filling consists of just 3 disks. 

Proof. Without loss of generality, we assume that s is directed vertically upward. We first assume that 
t is in the interior of D^(s). 

Consider the line st. Let us first assume that 7 intersects this line top-to-bottom (or tangentially) at t. 
In that case, we extend 7 forward by the semi-infinite ray starting at t, and backward by the boundary 
of Dl(s). We trace the resulting path backwards from infinity, stopping at the first intersection of the 
path with the part we have already seen, and thus forming a loop that does not contain t in its interior. 
We apply the Pestov-Ionin lemma to this loop, and find a unit-radius disk D contained in it. (See 
Figure W\) If D £ fil(P), then we are done. Otherwise, note that the loop does not cross the boundary 




Figure 7: (left) A unit disk D lies in the shaded area, (right) t is in lda(F). 

of P U Dl(s), so D £ P U Dl(s). Then we obtain a disk F £ lfil(s) by translating D upwards until it 
touches the forward chain FC(s), and we have t £ lda(F). 

Now we consider the case where 7 intersects the line st bottom-to-top at t. Let t' be the first point 
of intersection of 7 and the line st (along the path 7). If t' lies between s and t on the line st, then we 
can apply the argument above to conclude the existence of a disk F such that t' £ LDA(F) — but then 
also t £ lda(F). 

If, finally, t lies between s and t' on st, then we extend 7 by the semi-infinite ray starting in t. This 
ray must intersect the part of 7 from s to t' , and so again we have found a loop lying in P, fulfilling the 
requirements of the Pestov-Ionin lemma, and containing t. As above, there is then a disk F £ fil(P) 
with t £ lda(F). 

We now consider the case where t does not lie in the interior of Dl(s). Then t lies in a connected 
component P c of P \ Dl(s) different from LDA(s). Since -Dl(s) does not intersect the forward chain, P c 
lies entirely below Dl(s). Let t! denote the first point on 7 that is on the boundary of P c . We trace 
backward from t' a path along 7, and then along the lower semi-circle of Dl(s), until we reach a point 
that we have already seen. It forms a loop on which we apply the Pestov-Ionin lemma. Thus we find a 







disk D inside P U -Dl(s) that does not contain {t, t'}. If D £ FIL(P), then we choose F = D and we are 
done. Otherwise, we obtain F by translating D upward until it meets the forward chain, and we have 

t e lda(f). □ 



The lemma above does not give us a complete characterization of the reachable region, as we do 
not know yet whether we can reach the disks in fil(P) and lfil(P) tangentially. The following lemma 
addresses this issue. 

Lemma 8. Assume that the forward chain fc(s) does not intersect the interior of D^(s). (i) If F € 
fil(P), then there exists a counterclockwise configuration tangent to F that can be reached from s by a 
CS path, (ii) If F G lfil(P), then the first configuration on fc(s) tangent to F can be reached from s 
by a CSC path. 

Proof. The lemma is obvious when F touches the edge containing s, so we assume this is not the case. 
We first prove (i). If we draw a line segment upward from the leftmost point p of F until we meet the 
forward chain, we do not intersect Dl(s). It follows that p g lda(s). We consider a ray that starts at a 
configuration s' tangent to Dl(s). We start at s' = s and move s' counterclockwise along the boundary 
of Dl(s), so that the ray sweeps lda(s). Since p G lda(s), this ray must meet F at some point. When 
the ray first meets F, it is tangent to F, so we can reach the corresponding configuration by a CS path. 

We now prove (ii). We denote by r = (r, d) the first configuration on the forward chain that is tangent 
to F. Then r <G lda(s), so when we sweep the same ray as in the proof of (i), we meet F tangentially 
at some configuration r'. The arc between r' and r is inside P, so we have a CSC path starting from s 
and going through r' and r. □ 

We are now able to give the following characterization for the reachable region starting from a 
configuration on the side of P. It follows directly from Lemmas |6j[7j and [8] 

Proposition 9. Assume that s is a configuration on the boundary of P, oriented counterclockwise. Then 
any point in REACh(s) can be reached by a CSCS path. In addition, we have that: 

(i) If FC(s) intersects the interior of Dl(s), then REACh(s) = lda(s). 

(ii) 7/fc(s) does not intersects the interior of Dl(s), then REACh(s) = Ufgfil(P)ulfil(s) lda(P). 

In the characterization above, it seems that an infinite number of disks could possibly contribute 
to the boundary of the reachable region. In the following, we show that the contribution of the ldas 
along any edge can be reduced to at most two ldas. We focus on a particular edge /. Let d denote the 
counterclockwise direction along this edge. Then we order the counterclockwise configurations along / 
according to direction d, that is, for two such configurations Si = (si,d) and s 2 = (s 2 ,d), we say that 
Si =4 s 2 when Sist ■ d ^ 0. 

Lemma 10. Let Si — (si,d) and s 2 = (s 2 ,d) be two counterclockwise configurations on the same edge f 
of P, such that S\ =4 s 2 . Let hi = (/ii,d) denote the first counterclockwise configuration on f such that 
si hi and D L (hi) intersects FC(hi) \ {hi}. Then lda(s 2 ) C LDa(si) U LDA(hi). 

Proof. We first assume that fc(si) intersects the interior of -Dl(si), so that hi = si. Then lda(si) can 
be enlarged into a pocket, so by Lemma[3j we have lda(s 2 ) C lda(si). 

Otherwise, fc(si) does not intersect the interior of Pl(si). Thus, the disk Di(hi) touches FC(si). 
(See Figure|8]) If hi =^ s 2 , then LDA(hi) is a pocket, so by Lemma[3j we have lda(s 2 ) C LDA(hi). 

Finally, we assume that si =^ s 2 =<; hi. Let t denote a point in lda(s 2 ). If t is reached after an 
arc of Pz,(s 2 ) with length less than it followed by a line segment r 2 £, then it is clearly in LDA(si). (See 
Figure [8j left). On the other hand, if t is reached by an arc of L>l(s 2 ) with length at least n, followed 
by a segment r 2 i, we claim that t e LDA(hi). Let r be the point of P^(hi) such that the line rt is 
tangent to D^(hi) from the left. (See Figure[8j right.) We have to argue that the arc of _D^(hi) from hi 
to r lies in P. This follows from the fact that P^(hi) is obtained by translating L>l(s 2 ) along / until 
it touches fc(si), that the arc of P_l(s 2 ) from s 2 to r is in P, and that the arc from hi to r 2 is shorter 
than the arc from s 2 to r. □ 

Now we can show how to construct the reachable region from a configuration on the boundary. 
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Figure 8: Proof of Lemma [TO] 



Proposition 11. Let s be a counterclockwise configuration on the boundary of an n-sided convex 
polygon P. Then we can compute in 0(n) time a set bfil(s) of 0(n) disks, such that REACh(s) = 

Ufgbfi^s) LDA(F). 

Proof. We use the characterization of REACH(s) from Proposition [9j If FC(s) intersects the interior of 
Dl(s), then we just set BFIl(s) = {Dl(s)}. So in the remainder of this proof, we assume that FC(s) 
does not intersect the interior of Dl(s). 

If FIL(P) 7^ 0, then we first construct the contribution of the disks in fil(P). By Observation [4j we 
only need to find the unit disks whose centers are the vertices of the convex hull of the centers of the 
disks in fil(P). These disks are tangent to at least two edges of P, so their centers lie on the medial 
axis [UEl]. of P. We compute this medial axis in 0(n) time using an algorithm by Aggarwal et al. [I], 
and then check each edge on the medial axis to obtain these disks in O(n) time. 

We now observe that if -Dl(s) G fil(P), then we can set bfil(s) = fil(s) and are done. Thus, in 
the remainder of this proof, we assume that -Dl( s ) ^ fil(P), and explain how to find the contribution 
of lfil(P). 

Let d denote the first point on dP, starting from s in counterclockwise direction, such that d € Dl(s). 
(See Figure J9^1eft) .) Let us call a candidate configuration a configuration s' on the counterclockwise 
boundary of P with the property that Di(s') C PUDi(s) and such that Dl(s') is either tangent to two 
edges of fc(s), or is tangent to one edge of fc(s) and contains the point d. 

Consider an arbitrary edge / of P. Let Si denote the first counterclockwise configuration on / such 
that Dl(si) G lfil(P). When hi is as in Lemma 10 the ldas of all the disks in LFIL(P) that are 
tangent to / are contained in lda(si) U LDA(hi). So we only need to find Si and hi to construct the 
contribution of / to reach(s). 

We observe that S! and hi are candidate configurations: in fact, S! is the first candidate configuration 
on /, while hi is the last. 





Figure 9: On the left, P is shaded. On the right, P' is shaded, and its medial axis is dashed. 
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It remains to explain how to compute the candidate configurations efficiently. Denote by P' a simple 
polygon obtained by replacing the subchain of dP that goes counterclockwise from d to s with two or three 
edges, such that Dl(s) c P' . (See Figure fright).) The disk Dl(s'), for a candidate configuration s', 
must lie in P' and must touch dP' in more than one point. It follows that we can find the candidate 
configurations by first computing the medial axis of P' in 0(n) time using the algorithm by Chin et 
al. [TO], and then checking all edges of the medial axis. □ 



Proposition 11 shows that the reachable region for a configuration on the boundary of P is delimited 



by 0(n) disks. This bound is tight, as shown by the example in Figure 10 where Q(n) disks of lfil(s) 
contribute to the boundary of the reachable region. 




Figure 10: Example where f2(n) disks of lfil(s) contribute to the boundary of the reachable region. 
Here the disks of lfil(s) are centered along the dashed circle. The reachable region is shaded. 



4 Special left-right and right-left paths suffice 

Let s be a starting configuration. A canonical RL-start from s is a path from s to a configuration r on 
the boundary of P that begins with a right-turning arc of unit radius and continues with a left-turning 
arc of unit radius ending at r (and tangent to the boundary of P there) (Figure 111. 




Figure 11: A canonical RL-start. 

Note that for each edge / of P and for a given s, there are at most two canonical RL-starts from s 
ending on /. A canonical LR-start is defined analogously: it begins with a left-turning arc and continues 
with a right-turning arc. 

In this section, we show that for determining the reachability by paths in a convex polygon, it suffices 
to consider paths of a fairly special form. Namely, we show that a point is reachable if and only if it is 
directly accessible, or it can be reached by a path that begins with a canonical start. 
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Dubins [T2] showed that the shortest path of bounded curvature between two configurations in the 
plane is of type CSC or CCC. In the latter case, the middle arc has length more than it. (See Figure 12 ) 
We call these paths Dubins paths. 






(i) (ii) (iii) 

Figure 12: Three types of Dubins paths. 



Proposition 12. Let P be a convex polygon, let s be a starting configuration in P, and let t G P be 
reachable from s by a bounded- curvature path. Then t lies in the directly accessible region DA(s), or it 
can be reached by a path of one of the following forms: a canonical RL-start followed by a left-turning 
path (starting on a side of P), or a canonical LR- start followed by a right-turning path (starting on a 
side of P). 

Proof. Jacobs and Canny |15) showed that, in a polygonal environment, the shortest path of bounded 
curvature between two configurations is a sequence of Dubins paths. The final configurations of these 
Dubins paths (except for the last one) all lie on the boundary of the polygonal environment. So, if we 
denote by 7 a shortest path from s to i, then 7 can be written as a sequence 7 = 7172 • ■ • 7 m of m Dubins 
paths. When m ^ 2, we know that the final configuration q of 71 lies on dP. 

We handle three cases separately, according to the type of 71 (see Figure [i"2|) . If -Dl(s) and Dr(s) 
are contained in P, then DA(s) = P, so from now on, we assume that -Dl(s) or Dr(s) crosses dP. 





(a) 



(b) 



Figure 13: Proof of Proposition 12 case (i) 



Case (i). We assume that 71 is of type CCC. We denote by Ci, C2, C3 the three circle arcs such that 
71 = C1C2C3, and recall that C2 has length larger than ir. If C2 touches dP we are done, so from now 
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on we assume that C2 does not touch dP. Let D be the disk supporting C2. Without loss of generality, 
we assume that C\ and C3 turn counterclockwise and C2 turns clockwise. (See Figure [i"3^a).) Let cj, 
denote the center of Dl(s). We denote by Dg the disk obtained by rotating D by an angle 9 around cl- 
We define 9* as the smallest 9 £ (0, 2ir) such that 9 ^ tt/3 and D g \D touches dP, or 9 > ir/3 and Dg 
touches dP. Since -Dl(s) U -Dr(s) ^ P, such a 6>* exists. There is a path from s to q consisting of an 
arc of -Dl(s), an arc of Dg*, a line segment, an arc of D and C3. (See Figure 13 b).) This means that t 
can be reached by a canonical LR-start and a right-turning path. 







' Dg 


' 1 








Y 






CL 







(a) 



(b) 



(c) 



Figure 14: Proof of Proposition 12 case (ii). 



Case (ii). We assume that 71 = C1SC2, where G\ is left-turning, 5 is a segment and C2 is right-turning. 



(See Figure [14]) 

Let us first assume that m = 1, that is that 7 = 71 = C1SC2 and that C2 has length less than 7r. 
In this case, t lies in DA(s). Indeed, if C\ has length larger than tt or if t lies to the left of the directed 
line (! defined by s, then t g lda(s). (See Figure [T4|[c)). If C\ has length less than tt and t lies to the 
right of £', then t £ RDA(s), since 7 cannot enter Dr(s) because C2 has length less than tt. 

It remains to consider the case where m > 2 or the length of C2 is at least tt. Let £> be the unit 
disk tangent to Dl(s) and lying to the right of the segment S, and denote by Dg the disk obtained 
by rotating D by an angle 9 £ (0, 2tt] counterclockwise around c^. (See Figure 14 'a).) Let Cg be the 
clockwise arc of Dg starting at -Dl(s) n Dg and going clockwise until the first intersection point with 
the path SC2 , or returning to its starting point if there is no such intersection. Let 9* be the smallest 
value of 9 £ (0, 2tt] such that Cg touches dP. Since £>l(s) U Dr(s) <f_ P, such a 9* exists. Now q can 
be reached by a path consisting of an arc of Dl(s), an arc of Dg* , a line segment, and an portion of C%. 
(See Figure [HJb).) It follows that t can be reached by a canonical LR-start and a right-turning path. 






(a) 



(b) 



(c) 



Figure 15: Proof of Proposition 12 case (iii) 
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Case (iii). We assume that 71 = C1SC2, where C\ and C2 are right-turning and S is a segment. If 
the length of C2 is less than tt and to = 1 (see Figure [l5^ a)), then t £ Dr(s), and therefore t £ rda(s). 
We therefore assume that to ^ 2 or the length of C2 is at least tt . 

We denote by Dg the disk obtained by rotating D R (s) around the center of D L (s) by an angle 
9 e (0,7r/2]. (See Figure [l5|b).) When Dg intersects 71, we denote by Cg the arc of dDg that starts at 
Dl(s) n Dg and goes clockwise until it meets 71. Otherwise, we denote Cg = dDg. As before, let 9* be 
the smallest value of 9 £ (0, 2tt] such that Cg touches dP. Again, 9* exists since -Dl(s) U Dr(s) <f_ P. 
Now q can be reached by a path that consists of an arc of Dl(s), an arc of Dg*, a segment tangent to 
Dg* and 71, and a subpath of 71. (See Figure [il>|c).) This implies again that t can be reached by a 
canonical LR-start followed by a right-turning path. □ 

We can give a somewhat different characterization: 

Proposition 13. Let P be a convex polygon, let s be a starting configuration in P, and let t £ P be 
reachable from s by a bounded- curvature path. Then t is reachable by a path of the form CCSC'S. More 
precisely, t is reachable by a path of the form CS, or it is reachable by a path of the form CCSCS, where 
the two final disks touch the boundary of P, and the path goes through these touching points. 



5 Putting everything together 

In this section, we show how to construct the reachable region when s is an arbitrary configuration in P. 
We obtain it by combining the results in sections [3] and [4] We will prove the following: 

Theorem 14. Let P be an n-sided, convex polygon, and let s be a configuration inside P. Then the 
reachable region reach(s) from s inside P is delimited by 0{n) arcs of unit circles, and we can compute 
reach(s) in 0{n 2 ) time. 



Proof. Let t be a point in reach (s). By Proposition 12 either t is in da(s), or it can be reached 
after a canonical start. We only consider canonical RL-starts; the case of LR-starts can be handled 
symmetrically. 

The directly accessible region DA(s) is delimited by two circle arcs, which can be computed in 0{n) 
time by brute force. We determine the at most 2n canonical RL-starts by brute force, in 0(n 2 ) time. 
For each canonical start, by Proposition |11| we compute in Oin) time a set of 0(n) configurations on 
the side of P such that the union of their ldas form the reachable region after this canonical start. 

We have thus obtained a set of 0(n 2 ) configuration on the boundary of P such that the union of 
their ldas with da(s) is REACH(s). By Lemma [ToJ we only need to keep two such configurations per 
edge: the first and the last one. As we have only 0(n) arcs to consider, we can construct REACh(s) by 
inserting these arcs one by one, and updating the reachable region by brute force. As these arcs are arcs 
of unit circles, each one of them appears only once along the boundary of reach(s). So overall, it takes 
0(n 2 ) time. □ 
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